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Abstract A strong anisotropy in the longitudinal resistiv- 
ity of a 2D electron system has been observed at half-filled 
high Landau levels. We report on detailed Hartree-Fock cal- 
culations of the unidirectional charge density wave (UCDW) 
orientation energy induced by a tilted magnetic field. We find 
that stripes can orient both parallel or perpendicular to the 
in-plane field depending on the sample parameters and field 
strength. The close agreement between complex experimen- 
tal data on different sample geometries and our theoretical 
results strongly support the UCDW picture as the origin of 
the observed anisotropics in high Landau levels. 

1 Introduction 

Recently, a strongly anisotropic transport has been ob- 
served in high quality GaAs/Al2:Gai_2;As single hetero- 
junctions |l|,|,|,|,| at filling factors v = 9/2, 11/2, etc., 
and in 2D hole systems Q starting at v = 5/2. The 
origin of the magnetotransport anisotropy has not been 
firmly established yet. The most appealing interpreta- 
tion suggests that the 2D electron gas spontaneously 
breaks the translational symmetry by forming a unidi- 
rectional charge density wave, as predicted by Hartree- 
Fock theory Because of uncertainty about the relia- 
bility of this Hartree-Fock prediction, a special emphasis 
is placed on tests of its ability to explain experimental 
results on stripe orientation in tilted magnetic fields. 

Many-body RPA/Hartree-Fock calculations using a 
model, parabolic quantum well system showed that 
the stripe orientation relative to the in-plane field is not 
universal. In narrow quantum wells, stripes are expected 
to orient perpendicular to the in-plane field. For wider 
quantum wells, with two occupied subbands at zero tilt 
angle, a more complex behavior is predicted. Isotropic 
2D electron systems with half-filled lowest Landau level 
(LL) of the second subband can undergo a transition 
into a UCDW state due to level (anti)crossings induced 
by the in-plane magnetic field. The orientation of these 
stripe states can change from parallel to perpendicular 
as the strength of the in-plane field increases. 

In this paper we present calculations of the field-tilt 
UCDW anisotropy energy that include a self-consistent 
local-spin-density-approximation (LSDA) description of 
one-particle states in experimental sample geometries 
1^. Our results are compared to magnetotransport ex- 
periments on GaAs/Ala;Gai_a;As single hetero junctions 
with different 2D electron gas densities [||,|l§ and on a 
square quantum well sample with two occupied electric 
subbands [nil. 



2 Theoretical model and results 

The energy per electron of the UCDW state at fractional 
filling V* of the valence LL is given by 
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where a is the period of the UCDW state, e is the di- 
rection of charge variation, and An is the Fourier trans- 
form of the the guiding center occupation function at 
wave vector ril-K ja. In Hartree-Fock theory, the UCDW 
state energy depends only on a and e and the optimal 
UCDW is obtained by minimizing Eq. (|]) with respect to 
these parameters. In Eq.(|l|), t/(q) can be separated into 
direct, H{q), and exchange, -'^(q), contributions with 
1 
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where T^(q) is the RPA-screened effective Coulomb in- 
teraction 

To evaluate I^(q) we choose the in-plane component 
Bip of the magnetic field to be in the x-direction and 
use the following Landau gauge for the vector potential, 
A = (0, B±x — BipZ, 0). The one-particle orbitals are cal- 
culated using the self-consistent LSDA method and for 
any sample geometry, with the z-dependent confining 
potential, can be written as e^^^ / yjL^ipi^a{x — ("^kjz), 
where k is the wave vector which labels states within 
LL i, cr is the spin index, and i'^ — hc/eB±. The trans- 
lational symmetry responsible for LL degeneracy leads 
to a 2D wavefunction ipi^^(x,z) which is independent 
of the state label k, except for the rigid shift by i'^k 
along X-axis. This in turn leads to two-particle matrix 
elements of the Coulomb interactions with a dependence 
on state labels which is identical to that for the lowest 
LL of a zero-thickness 2D electron system provided the 
2D Coulomb interaction is replaced by the effective in- 
teraction V^(q). 

For an infinitely narrow electron layer the orbital 
effects of Bip vanish and the 2D Coulomb interaction, 
U(q), reduces toe-«'^'/2/g {LN{q^e'^ /2))^ 2Tre^£/e where 
L N (x) is the Laguerre polynomial and N is the orbit ra- 
dius quantum number. Starting from = 1, zeros of 
-^Af(9^^^/2) occur at finite q — q* , producing a zero in 
the repulsive Hartree interaction at wave vectors where 
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the attractive exchange interaction is strong. For the 
half-fiUed valence LL the corresponding UCDW state 
consists of ahernating occupied and empty stripes of 
electron guiding center states with a modulation period 
a w 27r/g*. 

In finite-thickness 2D systems the dependence of the 
effective interaction on wavevector magnitude q and ori- 
entation (j) relative to the in-plane field direction can be 
accurately approximated ||] by ^(q) — Vo{q) + V2{q) 
cos(20). At weak Bip, the isotropic term Va{q) has a 
wavevector-dependence similar to that of the effective 
interaction in the infinitely narrow 2D layer. The corre- 
sponding curve for the valence LL at v — 9/2 for a single 
heterojunction with 2D density rie = 2.67 x 10^"'^ cm~^ 
1^ is shown in the lower inset of Figure 1 . 

The non-zero anisotropy coefficient V2 (q) of the effec- 
tive interaction at Bip > 0, also shown in the inset, is re- 
sponsible for the non-zero UCDW anisotropy energy Ea, 
defined Q as the total Hartree-Fock energy of stripes ori- 
ented parallel with Bip minus the total energy of stripes 
perpendicular to Bip. The main plot of Figure 1 shows 
Ea calculated for single heterojunction samples. The two 
curves correspond to high and low density 2D electron 
systems with one occupied electric subband in perpen- 
dicular magnetic field. The subband wavefunctions are 
plotted in the upper inset. We found that stripes are ori- 
ented perpendicular to the in-plane field for the consid- 



ered interval of densities 



1.5-2.67X 10' Vm"'', con- 



sistent with available experimental data on single het- 
erojunctions |^,^,^. Similar agreement between theory 
and experiment is also abtained for other half-filled high 
LLs. 

In Figure 2 we present results for a 350 A wide GaAs 
quantum well with Ue = 4.6 x 10" cm"^ Q. At = 9/2 
and 11/2, and Bip — 0, the valence LL is the = LL 
of the second subband. The effective interaction Vo{q) 
has no zeros at finite wavevector, as seen in the inset 
of Figure 2, and the UCDW state is not expected to 
form, consistent with experiment [ pT| . An in-plane field 
Bip ^ 3 T induces a level anticrossing accompanied by 
a development of a minimum in Vo{q) at finite q. As 
discussed above, it is the minimum of the interaction 
energy at finite wavevector that opens the possibility for 
the formation of the UCDW state. The theoretical and 



experimental 1 1 1 in-plane fields corresponding to the on- 
set of the UCDW and transport anisotropy, respectively, 
are remarkably close. 

The main graph in Figure 2 indicates that, in this 



sample, the sign of Ea changes at B^. 



10 T. The 



stripes orient parallel with Bip at low fields while per- 
pendicular orientation is obtained in the high in-plane 
field region, consistent with the experimentally observed 
interchange of easy and hard current axes at similar Bip. 

We conclude that the close agreement between com- 
plex experimental data and theoretical results for all 
studied samples leave little doubt as to the UCDW origin 
of the observed transport anisotropies. 
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Fig. 1 Main graph; UCDW anisotropy energy for single 
heterojunction with 2D density Ue = 2.67 (solid line) and 
1.5 (dashed line) xlO^' cm~^. Upper inset: square of the 
normalized subband wavefunction for the two densities at 
Bip = 0. Lower inset: n = (thick line) and n — 2 (thin line) 
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Fig. 2 Main graph: UCDW anisotropy energy for a square 
quantum well, Tie = 4.6 x 10^' cm~^ and v = 9/2 (solid 
line) and v — 11/2 (dashed line). Inset: isotropic coefficient 
of the effective Coulomb interaction at = 9/2 and different 
in-plane magnetic fields. 

References 

1. H.L. Stormer et al, Bull. Amer. Phys. Soc. 38, 235 
(1993). 

2. M. P. Lilly et al, Phys. Rev. Lett. 82, 394 (1999). 

3. R.R. Du et al. Solid State Commun. 109, 389 (1999). 

4. W. Pan et al, Phys. Rev. Lett. 83, 820 (1999). 

5. M.P. Lilly et al, Phys. Rev. Lett. 83, 824 (1999). 

6. M. Shayegan et al, Physica E, 6, 40 (2000). 

7. A. A. Koulakov et al, Phys. Rev. Lett. 76, 499 (1996). 

8. R. Moessner and J. T. Chalker, Phys. Rev. B 54, 5006 
(1996). 

9. T. Jungwirth et al, Phys. Rev. B 60, 15574 (1999). 

10. M.P. Lilly and J. P. Eisenstein, unpublished. 

11. W. Pan et al, Phys. Rev. Lett, in press (2000). 



